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We carry out a systematic analysis of a pair of coupled qubits, each of which is subject to its own 
dissipative environment, and argue that a combination of the inter-qubit couplings which provides for 
the lowest possible decoherence rates corresponds to the incidence of a double spectral degeneracy 
in the two-qubit system. We support this general argument by the results of of an evolutionary 
genetic algorithm which can also be used for optimizing time-dependent processes (gates) and their 
sequences that implement various quantum computing protocols. 

The rapidly advancing field of quantum computing continues to bring about a deeper understanding of the basic 
quantum physics alongside prospects of the ground-breaking potential technological applications. However, one of the 
main obstacles on the way of implementing quantum protocols has so far been and still remains a virtually unavoidable 
environmentally induced decoherence. 

The list of the previously discussed decoherence suppression/avoidance techniques includes such proposals as error 
correction [1], encoding into decoherence-free subspaces (DFS) [2], dynamical decoupling/recoupling schemes (e.g., 
"bang-bang" pulse sequences) [3], and the use of the quantum Zeno effect [4]. 

However, none of the above recipes are entirely universal and their actual implementation may be hindered by 
such factors as a significant encoding overhead that puts a strain on the (still scarcely available) quantum computing 
resources (error correction), a stringent requirement of a completely symmetrical qubits' coupling to the dissipa- 
tive environment (decoherence-free subspaces), a need for the frequency of the control pulses to be well in excess of 
the environment's bandwidth (dynamical decoupling/recoupling) or an ability to perform a continuous high preci- 
sion measurement (quantum Zeno effect). While being more feasible in some (most notably, liquid-state NMR and 
trapped-ion) as compared to other quantum computing designs, the implementation of the above approaches might 
be particularly challenging in solid-state architectures. It is for this reason that augmenting the above techniques with 
a many-body physics-conscious engineering of robust multi-qubit systems and a systematic approach to choosing the 
optimal (coherence-wise) values of their microscopic parameters appears highly desirable. 

To this end, in the present work we explore yet another possibility of thwarting decoherence by virtue of permanent 
(albeit tunable) inter-qubit couplings. We find that, despite its being often thought of as a nuisance to be rid of, 
the properly tuned qubits' "cross-talk" may indeed provide for an additional layer of protection against decoherence. 
Specifically, we focus on the problem of preserving an arbitrary initial state ("quantum memory") of a basic two-qubit 
register during its idling period between consecutive gate operations. Moreover, the optimization method employed 
in our work can be further extended to the case of time-dependent two-qubit gates which can alone suffice to perform 
universal quantum operations. 

The previous analyses of the problem in question [5,6] have been largely limited to the symmetries of the underlying 
Hamiltonians and the parameter values corresponding to the presently available experimental setups [7-9]. As a result, 
they have not systematically addressed the issue of a possible role of the inter-qubit couplings in reducing decoherence. 
In contrast, our discussion pertains to a more general two-qubit Hamiltonian 

H= Yl [T.^aiB: + K{t)) + Jabal<E>a',], (1) 

a—x,y,z z— 1,2 

where each of the two qubits described by an independent triplet of the Pauli matrices ct^ (i = 1, 2) is subject to a 
local magnetic field comprised of the constant = (Ai,0, e.^) and fluctuating h^{t) — {0,0, hi{t)) components, the 
latter representing two uncorrelated (< hi(t)h2(0) >— 0) dissipative reservoirs. It is well known that the case of 
independent reservoirs present a significantly greater challenge than the greatly simplifying assumption of the highly 
correlated ones where the standard DFS can be readily found [2]. 

As regards the symmetry of the qubits' interaction, we restrict ourselves to the diagonal terms ct^ (8) errand the 
associated parameters Jab = Ja^ab- Thus, we exclude all the non-diagonal couplings between a\ and with a ^ b 
which, albeit possible in principle, do not normally occur in any of the known qubit designs. 

As we argue below, the preferred symmetries of the inter-qubit couplings can be found solely from the spectral 
analysis of the noiseless part of the Hamiltonian (1). This suggests that our results should be largely insensitive to 
the approximation used for treating the effect of the noisy reservoirs. 

To that end, we resort to the standard Bloch-Redfield (BR), i.e. a weak-coupling and Markovian approximation, as 
does most of the previous work on the subject [5]. Although the BR approximation is known to become potentially 
unreliable in the (arguably, most important for quantum computing-related applications) short-time (as compared to 
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the pertinent decoherence times) limit, the main advantage of the BR framework is a relative physical transparency 
of its results. Moreover, in light of the usual robustness of any symmetry-related conclusions, such as those to be 
described below, we anticipate that a more sophisticated analysis (e.g., akin that of Refs. [10]) would fully corroborate 
the BR results. 

In the standard singlet/triplet basis formed by the states: | tT>= (1,0,0,0)^, (| ti> +| iT>)/V2 = (0,1,0,0)^, 
I J.J.>= (0, 0, 1, 0)^, and(| ti> -| I1>)/V2 = (0, 0, 0, 1)"^, the noiseless part of the Hamiltonian (1) takes the form 
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where A = (Ai + A2)/V2, A" = (Ai - A2)/V2, e = ei + e2 and e" = ei - €2- 

The BR equations for the reduced two-qubit (size 4x4) density matrix p{t) takes a particularly simple form in the 
basis of eigenstates of Eq.(2) [5,11] 

Pnm{t) = -iiOnm Pnm{t) - ^ Rnmkl Pkl{t) (3) 

k,l 

where Wnm = [En — Em)/^ are the transition frequencies, and the partial decoherence rates 

Rnmkl = ^Im ^nrrk + ^nk ^Irrm ~ ^Imnk — ^knml (4) 

r r 

are given by combinations of the matrix elements of the relaxation tensor [11] 

Mmnk = -^SiuJnk) [(^l,lm(^l,nk + ^l,lm^l,nk] + JZ^ / S{w)(Lj " „ [o'z.ImO'z.nfc + O'z.JmO'z.nfc] (5) 

» 47r Jo w 

determined by the products between the matrix elements tr* computed in the eigenbasis of the noiseless Hamiltonian 
(2) and the spectral density of the reservoirs S{lo) — dte^'^*' < {hi{t), hi{0)} >{i = 1, 2). 

As in Refs. [5], we choose S{uj) to be Ohmic, S{uj) = acjcoth^, thus justifying the applicability of the Fermi's 
Golden rule-based expression (5) (hereafter we use the units ?i = fcs = 1). 

Obviously, the optimal choice of the Hamiltonian (2) would be the one that provides for the lowest decoherence 
rates composed of the matrix elements (5). The analysis of the expressions (5) reveals that their real parts tend to 
decrease monotonically with the increasing length J — \ J\ of the vector J ~ (J^, Jy, Jz)- However, in all the practically 
important cases (especially, in the before mentioned Josephson junction-based designs [7-9]) where the qubits emerge 
as effective (as opposed to genuine) two-level systems, an unlimited increase of J is not possible without leaving the 
designated "qubit subspace" of the full Hilbcrt space of the system. Besides, the coupling strength errors that are 
likely to occur in any realistic setup tend to increase with J as well. Therefore, in what follows we fix the length of 
the vector J and search for an optimal configuration of its components, allowing for either sign of the latter. 

Being proportional to the noise spectral density evaluated at frequencies corresponding to the transitions between 
different pairs of energy levels, the matrix elements (5), just like S{uj„i„) itself, attain their minimum values (which 
are proportional to the reservoirs' temperature T) at cJmn = 0. Thus, one might expect that the relaxation processes 
will be quenched and the decoherence rates will become suppressed, should some of the transition frequencies happen 
to vanish. 

Such a behavior would occur if any two of the energy levels of the Hamiltonian (2) became degenerate or, better yet, 
all the four levels of (2) became doubly (or even four-fold) degenerate. In this case, the lower pair of the degenerate 
levels could be viewed as an effective (" logical" ) qubit which is protected from decoherence by an energy gap separating 
it from the upper pair of levels, resulting in an exponential suppression of the relaxation rates at low temperatures. 
Then, by encoding quantum information into this logical qubit's subspace, one can greatly suppress its decay, any 
residual decoherence being solely due to pure dephasing controled by 5(0) ~ aT. 

Next, we apply this general argument to the practically important case of identical qubits (Ai = A2 = A/2) which 
are both tuned to the "co-resonance" point (ei = 62 = 0). 

We note, in passing, that a special significance of this parameter regime as a potentially most coherence-friendly 
one is also evidenced by the experiments on the Josephson charge-phase qubits [9] . 

At the CO- resonance point, the spectrum of Eq.(2) and the corresponding eigenvectors are given by the fol- 
lowing simpleexpressions: Ei = Jx + K,E2 = —{Jx + Jy + Jz),E3 = —Jx + Jy + Jz,E4 = Jx — K and 
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IXi >= {lAJy-Jz + K)/l^,l,Q)/^ 2K{K + Jy- J,), 1x2 >= (0,0 ,0,1),|X3 >= (-1, 0, 1, 0)/V2, \xa >= ih{Jy- 
- K)/A, 1,Q)/^J2K{K - Jy + J^, where K = {Jy - J^y + 2A^ , respectively. 

The only attainable incidence of double spectral degeneracy {Ei = £3, £'2 = E4 or Ei = E2,Es = E4) can occur for 
Jy.z > or Jy,z < 0, respectively. Resolving the above conditions, we find that the conjectured optimal configuration 
of the inter-qubit couplings must obey the equations 

J7* = 0, 2Jl^^rf = A^ (6) 

By taking into account the normalization of the vector J and Eq.(6) we finally obtain J"*"* = ±i(\/j^"+A2 — 
\J J"^ — A^), J°P* = ±i(\/ J2 + + \J J"^ — A^). Notably, this result shows that the conjectured optimal regime can 
only be achieved for sufficiently strong couplings (J > A). 

Having elucidated the physical content of the expected J°^*, we now confirm our predictions by applying a direct 
optimization procedure based on a genetic algorithm [12]. The latter starts out with random sets of parameters J 
which, in the language of Ref. [12], constitute an initial "population". In the course of the optimization procedure, 
the worst configurations are discarded and certain recombination procedures ("mutations") are performed on the rest 
of the population. The new solutions (" offsprings" ) arc selected according to the values of the chosen ffincss function. 
The iterations continue until the set of parameters converges to a stable solution representing the sought-after optimal 
configuration J°p*. 

As a fitness function we choose such a customary quantifier of the register's performance as purity P(t) = 
averaged over the solutions p^it) of Eq.(3) with the initial conditions /O-'(O) = |5'f„ >< 
given by the 16 product states %^ = >i ®\^b >2,{a,b = 1,...,4): |*i >= | i>,|4'2 >= | T>, 

1*3 >= l/\/2(| i> +1 T>),l*4 >= l/\/2(| i> +i\ T>). 

In Fig.l we plot the purity decay rate \dP{t)/dt\ as a function of Jy and for ei = Jx = and a = 10~^. For 
different values of J, the minima of of this two-parameter function collapse onto a hyperbola in the Jy — plane 
which appears to be described by Eq.(6) within the accuracy of our numerical solution of the BR equations (3). 

By varying the temperature, we observed that, consistent with our preliminary insight into the mechanism of 
the decoherence suppression, the effect of the interaction-induced double spectral degeneracy appears to be most 
pronounced at T <i; A, while at higher T the landscape shown in Fig.l flattens out, thus making the optimization 
less effective. 

We also observed that the use of an alternate performance measure, fidelity F = 3;^ T'r{p^(t)p^(t)} where 

PQ{t) represents the unitary evolution of an initial state p^{Q), yields the same results. 

In order to further illustrate our findings, in Fig. 2 wc present the behavior of a typical matrix element -R4123 as a 
function of a single parameter w introduced by the relations: = 0,Jy = 0.25Aw, = 1.98Au;. It can be readily 

seen from Fig. 2 that at the double degeneracy point w = 1 which corresponds to J = ^ Jj + Jy + J"^ = 2A the 

matrix element i?4i23 drops to its lowest value. 

By contrast, in the non-degenerate case exemplified by the Ising configuration Jx ~ Jy — 0, Jz = 2Aw the above 
dramatic dropis now absent, regardless of the value oi J = Jz- In Fig. 3 we contrast the purity P{t) computed for 
the optimal configuration with the results obtained in the Ising-, XY-, and Heisenberg- symmetrical, as well as non- 
interacting cases. These results suggest a systematic way of constructing a class of highly robust input states that 
undergo a slow decay governed by pure dephasing (no relaxation), as compared to the case of general encoding. 

For Jy^z > a natural candidate for an orthogonal pair of such states is presented by the (anti)symmetrical 
combination [ >= (|x2 > ±|X4 >)/V2 of the lower pair of the degenerate {E2 = E4) energy levels. The 

corresponding initial density matrix is characterized by the only non-zero entries p22{0) = ±/O24(0) = ±^42(0) = 
P44(0) = 1/2. 

In Fig. 4 we plot the decay rate of the purity dP{t)/dt as a function of the parameter w for the initial state |^^'''''« >. 
Note, that at the double degeneracy point {w = 1) the purity decay rate is reduced by several orders of magnitude. 

The high robustness of a degenerate ground state facilitated by the presence of the energy gap in our two-qubit 
system is somewhat reminiscent of the notion of " supcrcoherent" subspaces introduced in the abstract setting in Ref. 
[13]. However, despite a certain resemblance between the two, our implementation of a robust logical qubit is different 
in a number of aspects. 

For one, the construction proposed in Ref. [13] requires at least four physical qubits which interact by virtue of 
the Hamiltonian whose spectrum can be classified by the eigenvalues of the total angular momentum of the system. 
Then the degenerate singlet ground state behaves as an "error detecting code" that can only loose its coherence by 
absorbing energy from the reservoirs, the rate of which process appears to be exponentially suppressed for T — > 0. 
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In contrast, our logical qubit consists of only two physical qubits interacting via a generic (not necessarily spin- 
rotationally invariant) Hamiltonian, and its decoherence rate is limited by pure dephasing whose contribution remains 
linear in T even at the lowest temperatures. However, albeit not providing an equally strong protection for the logical 
qubit itself, in our case the double degeneracy of the two-qubit spectrum docs improve the gate performance for 
general encoding, as manifested by the gate characteristics averaged over different initial states. 

One can also identify the initial states whose decay can only get worse upon increasing the qubit coupling. These 
fragile states belong to the subspacc spanned by the complementary pair of cigcnstatcs \xi >, \x3 >■ 

Lastly, a word of caution is in order. Conceivably, any incidence of degeneracy may invalidate the use of the Golden 
rule-based Eq.(5) whose (sufficient) conditions of applicability requires that max\Rmnki\ ^ min\u}mn\- However, we 
found that our numerical solution remains stable even in a close vicinity of the double degeneracy point. 

The optimal configuration J°p* shows no abrupt changes upon varying the strength of the coupling a from its 
lowest values that do comply with the above sufficient conditions all the way down to the parameter range where 
minluJrnnl < max\Rrnnki\ < max\uJmn\- 

To summarize, in the present work we demonstrated that the properly chosen permanent inter-qubit couplings can 
provide a new means of protecting quantum registers against decoherence. While being relatively unexplored in the 
mainstream quantum information proposals, a similar idea has been in the focus of a number of scenarios exploring 
the possibility of a natural emergenceof logical qubits and the conditions facilitating fault-tolerant computations in 
strongly correlated spin systems. 

In the existing proposals [14] the prototype logical qubits are envisioned as topological ground and/or quasi-particle 
bulk/edge states of rather exotic chiral spin liquids. Thus, despite their enjoying an exceptionally high degree of coher- 
ence, the(topo)logical qubits require an enormous overhead in encoding, since creating only a handful of such qubits 
takes a macroscopically large number of interacting physical ones. Besides, the nearly perfect isolation of (topo)logical 
qubits from the environment can also make initialization of and read-out from such qubits rather challenging. Nonethe- 
less, our results show that a somewhat more modest idea of augmenting the other decoherence-suppression techniques 
with appropriate tailoring of the inter-qubit couplings might still result in a substantial improvement of the quantum 
register's performance. 

The numerical optimization method employed in this work can be further extended to the case of time-dependent 
two-qubit gates as well as beyond the Bloch-Redfield approximation and/or the assumption of the Ohmic dissipative 
environments. 

As regards the gate optimization, we find that in many cases the best performance of a two-qubit register can be 
achieved with a rectangular pulse Ja(t) = Ja6(T — t)Q{t) where the choice of the coupling parameters Ja facilitates 
the degeneracy of the lowest pair of energy levels. These results will be presented elsewhere [15]. 

The rapid pace of the technological progress in solid-state quantum computing gives one a hope that the specific 
prescriptions towards building robust qubits and their assemblies discussed in this work can be implemented in future 
devices. In this regard, very promising appear to be the phase [7], charge [8] and charge-phase [9] superconducting 
qubit architectures where, in principle, one can achieve any desired symmetry of the interaction terms in Eq.(l) by 
merely tuning the capacitive and inductive couplings between different Josephson junctions which implement the 
physical qubits. 

The authors acknowledge valuable discussions with J.E. Mooij and M. Dcvoret. 
This research was supported by ARO under Contract DAAD19-02- 1-0049. 
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FIG. 2. Partial dephasing rates exemplified by a typical matrix element i?4i23 as a function of the parameter w computed 
for T = 10"^ A, J = (0, 0.25Aw, 1.98Aw) (solid line) and J = (0, 0, 2Aw) (dashed line). 




FIG. 3. Purity as a function of time for zero coupling (J = 0, thin solid line); Ising (J = (J, 0, 0), dotted line and J = (0, J, 0), 
dashed line), XY (J = ( J/a/2, J/V^, 0), dash-double dotted line), Heisenberg (J — {J/VS, J/\/3, J/Vs), dash-dotted line), and 
the optimal configuration(J = J""*, sohd hne). Here J = 2A and T = 10~'^A 
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FIG. 4. Purity decay rate dP/dt as a function of the parameter w for J = (0, 0.25Aw, 1.98Aw)for T = 10 ^A (sohd line) 
and T = A (dashed line) . 
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